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شبیه‌سازی تغبیرات ریزساختار آلومینیوم ۱ در فرابند فشردن در کانال‌های هم‌مقطع زاویه‌دار 
با استفاده از روش مونت کارلو " 


3 ( 


7 هه ۱ ۰ ار بو و ۱3 
یکی از مهم‌ترین فرایندهای تغییرشکل پلاستیک شدیدء فشردن د رکانالهای هم‌قطع زاویه‌دار است. این روش دارای مزایای زیادی از جمله 
ریزدانگی بسیار شدید است. در این پژوهش, فرایند فشردن در کانال‌های هم‌طع زاویه‌دا رآلیا ژآلومینیوم 107۱ توسط نرم‌افزا رآباکوس 
شییه‌سازی و توزی عم کرنش در نمونه استخحراج شد. سپس ریزسانعتار نمونه آنیل شده پس از فشردن د رکانال‌های هم‌مقطع زاویه‌دار: توسط 
روش مونت کارلو» در محیط نرم‌افزار متلب» شبیه‌سازی و متوسط اندازه دانه نهایی پیش‌بین یگردید. نتایج شبیه‌سازی‌ها نشان می‌دهد که با 
افزايش ناهمکنی توزی مکرنش در سطح نمونه فشرده شده در کانال‌های هم‌قطع زاویه‌دار ناهمکنی اندازه دانه‌ها در ریزساعتار نهایی پس از 


آنیل نیز بیشتر حواهد شد. همچنین با افزايش زاویه کانال قالب, متوسط اندازه دانه نهایی پس ا زآنیل افزایش می‌یابد. 
واژه‌های کلیدی آلیاز آلومینیوم ۱ فرایند فشردن در کانال‌های هم‌مقطع زاویه‌دان اجزاء محدود. مونت کارلو اندازه دانه. 


۲ ۱۵۱۱۱6 امین ۱ ۸۵۵061 ۵ ممتاا۲۵ ۲ اما موم م۱۷۲۱ عط ۵۶ حمتامابصته 


27۱ ۱۷۲۵۳۸۸۵ وصزون۱ ۳۵۶910۵ 


م5206 ,۱۷۲ .ظ 060 ,۳۲ ٩.‏ ۸0۷/۵1061 ,۲۱ .1۷۲ 12۳67 .1۷۲ 


و( 
کمیوا0ع 0۳۳۵0۲]ع0 عاکهام 56۷۵۳۵ ۱۳۵۵0۳۱۵۳ ۵5۲ 66 0۳8۵ 15 (حلص) وصنعکع ۲ج موجه مه آمبو۲ 
۴ 2۲۵6656 ط 2۱6 ,ناه ۲6۹۵۱ 6 3 ۳۵2۵۵ 0و 6065۲ که مینک کعوهار ۳۵ ۱۳۱ 
۵۵ ۵۷۵۱۵۵۱۵0۰ کز (0باطزکن اک 6 ۵0۳0۵ 5۵/۵۵۲۵ دنام‌ضشظد وطاعب 0عامابونکی و ۸۸۵6061 
۵ و۱۱ 0عمآیعآه6 5 5۱26 9۲۵ 0۷۵۲۵96 626 0۳00 2۳۵06160 کا و0060 06 عا5۵ 2۲6۹560 126 ]0 
۱۱/۱۵0۳۱۵9۵ 5۳۵۲ ۱6 ۵5 ۱۵۲) عا۱۱06۵ کاآیوع۲ ۲۳۱6 50۵]۵۷۵۲۶۰ ۱۸۲ ]هن 0اه مه ۱ از مواه ماه 
660۳۵ ۱6 05 ,1۷0۲۵0۷6۲ ۳۵۵0۱6۵5۰ ۳۵۵0۱60 م0۵ 06 0 میک 9۳01۲ 16 ر5ع۱6۳۵۵5[ 

[۱۱6۲5600۳ 0۳۵16 1۳6۲۵۵565, 616 0۷6۲۵96 9۲۵1۳ 51276 15 0: 


۳۸۰ ۶۱26 صنفی زماتمن ماع1۷۲۵ :۳۵۸۳ :۸۵۸6061 16۷۵۲۵8 


نسخه نخست مقاله در تاریخ ٩۲/۱۲/۱۱‏ و نسخة پایانی آن در تاریخ ۹۵/۳/۲۲ به دفتر نشریه رسیده است. 
(۱) دانشجوی کارشناسی ارشد, دانشکده مهندسی مواد دانشگاه علم و صنعت ایران. 
(۲) نویسنده مسئول: استاده دانشکده مهندسی مواد دانشگاه علم و صنعت ایران. 22۲62۵)۵105].26.1۳ 
(۳) استاد , دانشکده مهندسی مواد. دانشگاه علم و صنعت ایران. 
(۶)استادیار: دانشکده مهندسی مواد. دانشگاه علم و صنعت ایران. 
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معد مه 


امروزه روش‌های تغییرشکل پلاستیک شدید با سرعت 
زیادی در حال گسترش بوده و بسیاری از آن‌ها در 
مرحله گذر از حالت آزمایشگاهی به صنعتی می‌باشند. 
این روش‌ها توانایی تولید آلیاژهای دارای ساختار بسیار 
ریز در مقیاس نانومتر را فراهم ساخته‌اند. یکی از مهم- 
ترین فرایندهای تغییرشکل پلاستیک شدید. فشردن در 
کانال‌های هم‌مقطع زاویه‌دار (۳۵۸۵۲ - اعحصصقط لقنوم 
0۳6558 ۲قآنا20۵) می‌باشد. از مهم‌ترین مزیت‌های 
فرایند ۳۵۳ ایجاد برش ساده تقریبا یکنواخت» شدید 
و زاویه‌دار درحجم بیلت می‌باشد. این مزیت یک نقش 
کلیدی را در تکامل ساختار برای کرنش‌های بزرگ ایفا 
می‌نماید. به این صورت که تنش‌های برشی بوجود آمده 
در ماده باعث بالا رفتن چگالی نابجایی‌ها در آن 
می‌شود. اين نابجایی‌ها. مرزدانه‌هایی با زاویه کوچک را 
تشکیل می‌دهند که در نهایت با تبدیل شدن به 
مرزدانه‌های با زاویه بزرگ ردیف‌هایی از دانه‌های ریز 
را بوجود می‌آورند. کوچک شدن اندازه دانه‌ها موجب 
بالا رفتن خواص مکانیکی و سوپرپلاستیکی مواد 
می‌شود. از دیگر مزایای مهم این روش تکرارپذیر بودن 
آن است که در این صورت امکان رسیدن به اندازه 
دانه‌های بسیار ریز میسر می‌گردد [1]. 

همان‌طور که گفته شد. در حین فرایند تغییرشکل» 
خواص مکانیکی مواد به صورت گسترده‌ای تغییر 
می‌کند. ماده تغییرشکل یافته از لحاظ ریزساختاری» ریز 
دانه شده و مقدار زیادی انرژی کرنشی در آن ذخیره 
می‌شود. در صورتیکه پس از تغییرشکل, ماده در دمای 
بالایی حرارت داده شود (آنیل شود). فرایندهای فعال 
شده حرارتی, مکانیزم‌هایی را فراهم می‌آورد که می‌تواند 
باعث حذف عیوب يا آرايش دیگری با انرژی کمتر 
شود. دو فرایند اصلی برای ترمیم ریزساختار وجود 
دارد: بازیابی و تبلور مجدد. اين فرایندها با دما فعال 
می‌شوند و در دو نوع استاتیکی و دینامیکی رخ می‌دهند 
[2]. 

فرایندهای بازیابی و تبلور مجدد بطور گسترده‌ای 


در آزمایشگاه‌های تحقیقاتی مورد بررسی و تحقیق قرار 


نشريه مهندسی متالویژی و مواد 


شبیه سازی تغییرات ریزساعتا رآلومینیوم... 


گرفته‌اند. اما امروزه گرایش به شبیه‌سازی‌های 
کامپیوتری برای پیش‌بینی تغییرات ریزساختاری با 
سر عفر تال افتوایش ابش کتاویزه یه ستازق 
کامپیوتری با استفاده از روش‌های عددی جهت مدل- 
جافی ود دنو فلور مه اناوت سا ۱/۸۰ 
شروع شد. روش‌های عددی رایج مورد استفاده عبارتند 
از؛ (۱) روش ورتکس. (۲) مدل زمینه فازی. (۲) روش 
اتوماتیک سلولی» (4) مدل مونت کارلو [3]. 

یکی از رایج‌ترین روش‌های عددی برای مدل 
کردن تبلور مجدد که در دو دهه اخیر مورد توجه قرار 
گرفته است. روش‌های مبتنی بر روش مونت کارلو 
شکور ال ۱۹/۳ او سم وا اسفاوه امن 
مونت کارلو در زمینه پیش‌بینی رشد دانه توسط 
اسرولویتز و همکارانش ارائه گردید [4]. 

از سال ۱۹۸۲ تاکنون تحقیقات زیادی توسط 
محققین مختلف جهت شبیه‌سازی تبلور مجدد و رشد 
دانه توسط مدل مونت کارلو صورت پذیرفته است -5] 
[18. تحقیقات انجام شده در راستای شبیه‌سازی دو و 
سه بعدی رشد دانه به صورت نرمال, غیر نرمال و رشد 
دانه در حضور رسوبات ریز می‌باشد. در مجموعه 
تحقیقاتی که تاکنون انجام شده سعی گردیده جهت 
تطبیق نتایج حاصل از مدل با شرایط واقعی رشد دانه 
این مدل اصلاح شود. نتایج بررسی‌ها نشان می‌دهد که 
سینتیک و نحوه هندسی رشد دانه در شبیه‌سازی رشد 
دائه با استفاده از این مدل با شرایط واقعی مطابقت دارد. 
در تحقیقات زیادی که بر روی پدیده تبلور مجدد انجام 
گرفته. فرض شده است که عامل محرک برای تبلور 
مجدد. یعنی انرژی ذخیره شده. ناشی از تغییرشکل 
همگن می‌باشد. به این منظور مقدار انرژی ذخیره شده 
ای از کر وا دک وی 
انرژی مرزدانه در ماده فرض گردیده است [17]. اما 
شرایط واقعی فرایند به همراه شکل هندسی قطعه کار 
اغلب منجر به توزیع ناهمگن انرژی ذخیره شده و در 
نتیجه منجر به جوانه‌زنی در مکان‌های خاص (توزیع 
غیر تصادفی) درون قطعه‌کار خواهد شد. به عنوان مثال 


سال بیست و هشتم شماره یک. ۱۳۹۵ 


منصوره طاهری - محمدرضا ابوطالبی - سید حسین سیدین- باقر محمدصادقی ۳۹ 


در مورد اکستروژن برای تولید پروفیل‌های آلومینیوم. 
بنابراین. تهیه الگوریتمی برای محاسبه توزیع اندازه دانه 
در آلیاژ تغییرشکل یافته ناهمگن و سپس تبلور مجدد 
یافته ضروری می‌باشد. در میان تحقیقات انجام شده 
جهت شبیه‌سازی ریزساختار تبلور مجدد يافته به وسیله 
مدل مونت کارلو فقط چند محقق انرژی ذخیره شده 
توسط نابیجاپی‌ها را در ماده بطور غیرهمگن در نظر 
گرفته اند [6,16,18]. 

تاکنون تحقیقات تتوری و آزمایشگاهی زیادی بر 
روی فرایند فشردن در کانال‌های هم‌مقطع زاویه‌دار و 
تاثیر پارامترهای مختلف بر آن انجام شده است. اماء به 
تکامل ریزساختاری بررسی‌های اندکی در مورد تاثیر 
این پارامترها بر پالایش ریزساختار صورت گرفته است. 


دلب پیچیدگی مکانیزم و پدیده‌های ناشناخته مربوط به 


بنابراین» در این‌باره به مطالعات بیشتری نیاز می‌باشد 
[19]. 

با توجه به مطالبی که ذکر شد. این پژوهش بادو 
هدف کلی صورت گرفته است: در مرحله اول فرایند 
فشردن در کانال‌های هم‌مقطع زاویه‌دار آلباژ آلومینیوم 
۱ توسط نرم‌اف زار اجزاء محدود آباکوس 
(۸3۸075) شبیه‌سازی شده و تاثیر برخی پارامترهای 
فرایند در آن مورد بررسی قرار گرفته است. در مرحله 
بعد. ریزساختار نمونه فشرده شده در کانال‌های 
هم‌مقطع زاویه‌دار پس از آنیل توسط روش مونت کارلو 
شبیه‌سازی شده و متوسط اندازه دانه نهایی پیش‌بینی 
می‌شود. به این منظور. از مرتبط کردن دو روش اجزاء 
محدود در محاسبه توزیع کرنش و مونت کارلو جهعت 
پیش‌بیشی انداوه فانه انتفاده شنده انس 


روش تحقیق 
در این بخش یک مدل و راه حل تشوری جامع ارائه 
خواهد شد که در آن پس از تعیین توزیع کرنش در 
نمونه فشرده شده در کانال‌های هم‌مقطع زاویه‌دار توسط 
شبیه‌سازی اجزاء محدود. توزیع انرژی ذخیره شده 
ناشی از تغییرشکل محاسبه شده و در ادامه ریزساختار و 


توزیع اندازه دانه در نمونه بعد از آنیل تعیین می‌شود. 
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این راه حل شامل کوپل روش اجزاء محدود و مونت 
کارلو می‌باشد. روندنمای کلی مدل ارائه شده در شکل 
)۱( نشان داده شده است. 


پیش بیلیریزساختارنهیی 
(توسط مدلمونت کار 


شک ۱ روتدنای کلی عللارافدشده 


شبیه سازی احزاء محدود فرایند فشردن در کانال‌های 
هم‌مفطع زاویه‌دار. شبیه‌سازی توسط ‏ نرم‌افزار 
ان‌نا:۸۸0/۹/۲ و به صورت سه بعدی انجام شد. 
در این راستا مجموعه سنبه و قالب به صورت صلب 
گسسته (3810 عاع56ز1) و بیلت به صورت انعطاف- 
پذیر (1۱6/0۳۳020[16) در نظر گرفته شد. 

قالب فرایند فشردن در کانال‌های هم‌مقطع 
زاویه‌دار با مقطم ۱۸*۱۶ میلیمتر مربع» زاویه کانال ٩۰‏ 
درجه و زاویه گوشه ۳۰ درجه و نیز نمونه‌ای با مقطع 
۶ میلیمتر مربع و طول 11 میلیمتر طراحی شد. 
شکل (۲ قالب‌های رسم شده در نرم‌افزار را نشان 
می‌دهد که به ترتیب از چپ به راست دارای زوایای 
کاثال ۵۰ ۱۲۳و ۱۳۵ وهی باشتا: 


شکل ۲ قالب‌های فرایند فشردن در کانال‌های هم‌مقطع زاویه‌دار 


طراحی شده در نرم‌افزار آباکوس 


در این شبیه‌سازی مقدار ضریب اصطکاک ۰/۰۹۵ در 
نظر گرفته شد. علاوه بر این سرعت حرکت سنبه یک 
میلیمتر بر ثائیه در نظر گرفته شد. 

در ادامه به منظور شبیه‌سازی فرایند فشردن در 
کانال‌های هم‌مقطع زاویه‌دار در پاس دوم با کپی کردن 
مدل پاس اول و جابجایی قالب و سنبه در ماژول 
7 و همچنین عوض کردن شرایط مرزی در 
ماژول 1.020 شرایط مسیر 3 پیاده سازی شد. لازم به 
ذکر است که در پاس دوم در ماژول 020 شرط اولیه 


محاسبه انرژی فخیره شده در تواحی مخاف تمونه 
فشرده شده در کانال‌های هم‌سفطع زاویه‌دار. همانطور 
که گفته شده تشییرشکل پلاستیکی شدیده یک ساغتار 
زیردانه با دانسیته نابجایی بالا در ماده بوجود می‌آورد 
ار مخت باه ات از کارت فک زاس 
حجم. 1 از رابطه زیر محاسبه می‌شود [19,20]: 
)0( 1-2 
که در این رابطه 1 متوسط اندازه زیردانه و 7 
انرژی بر واحد سطح یک مرز زیردانه است که در مرتبه 
۵ ول بر متر مربع می‌باشد. باید ذکر نمود که اندازه 
زیردانه با افزایش کرنش کاهش می‌یابد و گیل سویلانو 
و همکاران [21 رابطه زیر را با انجام تست بر روی 
فلزات مختلف بدست آورده‌اند: 
)۲( 35+172) <ظ 


آنجایی که در فرایند فشردن در کانال‌های هم‌مقطع 
زاویه‌دان توزیع کرنش در طول نمونه یکنواخت نمی- 
باشد (حداقل در پاس‌های اولیه» بنابراین انرژی ذخیره 
شاه نتفای با اند ات مس اش سای 
محاسبه انرژی در نقاط مختلف. سطح مقطع نمونه به 
۰ منطقه مساوی تقسیم و در هر منطقه متوسط مقدار 


کرنش اعمالی محاسبه شد. در واقع فرض شد که در هر 


نشریه مهندسی متالویژی و مواد 


شبیه سازی تغییرات ریزسااعتا رآلومینیوم... 


منطقه تغییرشکل همگن می‌باشد. در ادامه میزان انرژی 
ذخیره شده در هر منطقه توسط روابط مذکور محاسبه 
شد. قبل از انجام شبیه‌سازی تبلور مجدد توسط مدل 
مونت کارلو باید انرژی ذخیره شده را به همه مکان‌های 
شبکه ورودی اولیه به مدل اعمال نمود. 


پیش‌بینی ریزساختار نهایی توسط مدل مونت کارلو 
تبلور مجدد و رشد دانه. جهت شبیه‌سازی ریزساختار 
نمونه فشرده شده در کانال‌های هم‌مقطع زاویه‌دار پس از 
آنیل» از یک مدل مونت کارلو دو بعدی تبلور مجدد و 
رشد دانه استفاده شد. این مدل در محیط برنامه نویسی 
نرم‌افزار متلب (ظ۷۲۸۲1۸۳) پیاده‌سازی شد. در این 
راستاء یک شبکه دو بعدی مربعی با ابعاد ۱۰۰۱۰۰ 

سپس به هر مکان از شبکه یک عدد تصادفی 
[ی0) ب1]:٩‏ نسبت داده شد که نشان دهنده جهت هر دانه 
انجام شده این عدد قابلیت ایجاد رشد نرمال دانه را 
دارد[14,13,9]. شایان ذکر است مکان‌های دارای :5 
ایجاد ریزساختار اولیه مناسب قبل از تبلور مجدد. شبکه 
باید. علاوه بر انرژی ذخیره شده در بخش قبل. دارای 
اندازه اولیه مناسب نیز باشد. همچنین باید دانسبته 
محاسبه و در شبکه اعمال گردد. روش‌های مورد 
استفاده برای آماده‌سازی شبکه اولیه مناسب و همچنین 
الگوریتم مونت کارلو پیاده‌سازی شده در ادامه توصیح 


داده شده است. 


حول شک ارات با معراط اند از اه ما مه 
متوسط اندازه دانه آلیاژ آلومینیوم 1۰7۱ پس از انجام 
فرایند فشردن در کانال‌های هم‌مقطع زاویه‌دار در پاس 
اول در حدود ۱۸ میکرومتر و در پاس دوم حدود ۰/۵ 
میکرومتر می‌باشد [22]. برای ایجاد اندازه دانه مناسب 
می‌بایست از الگوریتم رشد نرمال دانه مونت کارلو 
استفاده نمود[13,14]. روندنمای الگوریتم رشد نرمال 
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دانه در روش مونت کارلو مورد استفاده در شکل (۳) 
نشان داده شده است. 


ت_ 


به ازای 
10*160 ۱21 


شکل ۳ روندنمای الگوریتم رشد نرمال دانه در روش مونت کارلو 


مورد استفاده 


هنگامی که ۱۰۰*۱۰۰ مرتبه مراحل فوق (به 
عبارت دیگر تغییر جهت) انجام پذیرد. یک گام مونت 
کارلو (۷805) انجام شده که واحد زمانی شبیه‌سازی 
می‌باشد. پس از اتمام هر گام شبکه حاصل مورد آننالیز 
قرار گرفته و متوسط اندازه دانه در آن محاسبه می‌شود. 
الگوریتم فوق تا زمان دستیابی به متوسط اندازه دانه 
مورد نیاز تکرار می‌شود. 
ریزساختار حاصل از آنیل می‌بایست توزیع انرژی 
ذخیره شده مورد محاسبه از روش مذکور به شبکه 
مونت کارلو منتقل گردد. بنابراین در این مرحله. شبکه 
مونت کارلو نیز به ۱۰ منطقه تقسیم شد که هر یک 
شامل ۱۰ ردیف ۱۰۰ مکانی یعنی ۱۰۰۰ مکان می‌باشد. 
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سپس به هر یک از مکان‌های یک منطقه. انرژی داخحلی 
مطابق با آنچه قبلاً ذکر شد محاسبه گردیده تخصیص 
داده می‌شود. از شبکه حاصل به عنوان ورودی الگوریتم 
تبلور مجدد مونت کارلو استفاده شد. 


مدل جوانه‌زنی. حال پس از انتقال توزیع انرژی ذخیره 
شده به شبکه مونت کارلو. شبیه‌سازی جوانه‌ها مطابق 
انرژی هر منطقه انجام می‌گیرد. از آنجایی که جوانه‌زنی 
با حرکت نابجایی‌ها» که فرایتدی فعال شونده با دسا 
است. روی می‌دهد. دانسیته جوانه‌زنی را می‌توان 
براساس رابطه‌ای بیان نمود که از یک عبارت برای نرخ 
ثابت فرایند. یک عبارت برای عامل محرک و یک 
ارت تفای مک ان ات رم اروت سک له 
است [16]: 


۳( | جک اعد رطع 


کهفر ان )ارت ناگی بط ان ری یه 
شده 1 انرژی ذخیره شده بحرانی که در مقادیر کمتر 
از آن جوانه‌زنی رخ نمی‌دهد» ,2 انرژی اکتیواسیون 
برای جوانه‌زنی و دمای آنیل می‌باشد. از آنجایی که 
در نمونه فشرده شده در کانال‌های هم‌مقطع زاویهدان 
انرژی ذخیره شده در ۱۰ منطقه محاسبه شده یکنواخت 
نبوده بنابراین دانسیته جوانه‌زنی نیز در هر ناحیه 
متفاوت خواهد بود. مقدار انرژی اکتیواسیون برای آلیاژ 
آلومینیوم ۲۰۱۱ برابر با ۱۵۰ کیلوژول بر مول می‌باشد 
[23]. 

شایان ذکر است که در اینجا از کاهش انترژق 
مربوط به بازیابی در ابتدای آنیل صرف نظر شده است. 
در واقع فرض شده است که در فرایند آنیل نرخ گرم 
شدن به اندازه‌ای زیاد بوده که فرصت بازیابی از فلز 
گرفته شده است بنابراین رابطه مربوط به دانسیته 


جوانه‌زنی به صورت غیر وابسته به زمان می‌باشد. 
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مدل مونت کارلو جهت شبیه‌سازی تبلور مجدد و 
رش دانه. روندنمای شبیه‌سازی تبلور مجدد با استفاده 


از روش مونت کارلو در شکل() نشان داده شده انبنتگا: 


داده‌های درودی 


آماده سازي شبکه اولیه 


محاسبه و اعمال جوانه ها 


پذیرتی غدم پذیرتی رند داند 


1۳0 1 


شکل ۶ روندنمای شبیه‌سازی تبلور مجدد با استفاده از روش 


مونت کارلو 


شبکه اولیه ورودی به حلقه اصلی برنامه شبکه‌ای 
با ابعاد ۱۰۰۱۰۰ با متوسط اندازه دانه۰/۵ میکرومتر 
می‌باشد که جهت گیری دانه‌ها عددی بین ۱ تا ۶۸ است. 
همچنین به هر مکان از شبکه. یک انرژی داخلی ناشی 
از تغییرشکل نیز اختصاص داده شده است. از آنجایی 
که متوسط اندازه دانه اولیه ۰/۵ میکرومتر می‌باشد. 
بنابراین در اینجا؛ اندازه هر مکان از شبکه برابر با ۰/۲۵ 
میکرومتر فرض شد. بنابراین با یک شبکه ۰۰۱۰۰ 
سطح 1۲۵ میکرومتر مربع را می‌توان شبیه‌سازی نمود. 

برای الگوریتم جوانه‌زنی از مدل سرعت ثابت 
استفاده شد. در این مدل. همانطور که قبلا توضیح داده 
شده تمامی جوانه‌ها در ابتدای آنیل اعمال می‌شود. در 
واقع برای بعضی از مواد. دانسیته جوانه‌زنی از ابتدای 
شروع آن ممکن است بسیار سریع کاهش یافته و به 
صفر برسد. به‌طوری که می‌توان فرض نمود تمام 
جوانه‌ها در لحظه شروع تبلور مجدد ایجاد شده و پس 
از آن هیچگونه جوانه‌زنی صورت نمی‌گیرد [2] 


نشریه مهندسی متالورژی و مواد 


شبیه سازی تغییرات ریزساختا رآلومینیوم... 


محاسبه دانسیته جوانه‌ها برای هر منطقه از شبکه 
مونت کارلو مطابق روابط ذکر شده می‌باشد. جهعت 
ایشا حغوانه‌هاه به ازاق هه خوال سک مکان از شیک 
به‌طور تصادفی انتخاب و یک عدد تصادفی در بازه 1٩‏ 
تا 1۶ به آن نسبت داده شد تا آن مکان بعنوان جوانه 
بدون انرژی داخلی درنظر گرفته شود و مقدار انرژی 
ناه ان کاخ یط بت ات هتفر وورزط, 
گرفته شود [9,24,25] 

پس از مرحله فوق, الگوریتم رشد دانه مونت 
کارلو به شبکه جوانهریزی شده اعمال گردید. 
روندنمای این الگوریتم همانند شکل (۳) می باشد. با 
این تفاوت که: 

۱- در محاسبه انرژی هر مکان علاوه بر در نظر گرفتن 
انرژی مرزدانه‌ها» انرژی ذخیره شده ناشی از 
تغییرشکل هر مکان نیز منظور گشته و از رابطه زیر 
محاسبه شد [9,16,18]: 

(( (5-,۵) ۷۶+ )1 - نا 

1 
در این رابطه ‏ انرژی مرزدانه است که بر واحد سطح 
گزارش می شود. (,5-,8)0 تابع پله‌ای است که در 
شرایط ,5 < ,4 برابر با یک و در ,5 > ,0 برابر با 
مق و 6 برای بتا 8۸ صری که لام ات 
دانه‌های اولیه تبلور مجدد نیافته می‌باشد. 1 مقدار 
ثابت مثبتی است که بیانگر انرژی ذخیره شده هر مکان 
از شبکه می‌باشد و مقدار آن از رابطه انرژی ذخیره شده 
مذکور و با توجه به مقیاس طولی شبکه محاسبه 
می‌شود[15,25]. 

۲- از تغییر جهت‌هایی که منجر به تبدیل دانه تبلور 
مجدد يافته به دانه تبلور مجدد نیافته می‌شود. 
جلوگیری به عمل آمد. به این منظور در صورتی 
که مکان انتخابی جز مکان‌های تبلور مجدد نیافته 
باشد (,0 > )۰ مجاز به تغییر جهت به کلیه 
هیا ها ی با شاوی تفن کته باب زا سره 
یکی از همسایه‌های تبلور مجدد يافته پذیرفته شود. 


رشد دانه تبلور مجدد یافته به درون دانه تبلور 
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مجدد نیافته رخ داده است (که با کاهش انرژی 
شیر ۵ نله هر اه امیبت )ه ق دی تالا کته تفر 
جهت به یکی از همسایه‌های تبلور مجدد نيافته 
پذیرفته شود رشد دانه عادی رخ داده است (که با 
کاهش انرژی مرزدانه همراه است). اما در صورتی 
که مکان انتخابی جز مکان‌های تبلور مجدد یافته 
باشد (,0 < )۰ تنها مجاز به تغییر جهت به یکی 
از همسایه‌های تبلور مجدد يافته می‌باشد. 


نتایج و تفسیر آن‌ها 

تتایج شبیه سازی‌های نمونه تحت فرایند فسردن در 
کانال‌های هم‌مقطع زاویه‌دار. اولین شبیه‌سازی فرایند 
فشردن در کانال‌های هم‌مقطع زاویه‌دار بر روی آلیاژ 
الومینیوم ۱۰۱۱ در قالبی با زاویه کانال ٩۰‏ درجه و 
زاویه گوشه ۳۰ درجه در دمای محیط و برای یک 
پاس؛ در نرم‌افزار آباکوس صورت پذیرفت (2561)). 
تاش بارآمت‌های اقعای تور تفن زر( ان واد: 
شده است. در فرایند فشردن در کانال‌های هم‌مقطع 
زاویه‌دار بدلیل تغییرشکل‌های زیاد و کرنش بالا در 
قالب ٩۰‏ درجه و همچنین بالا بودن اثر اصطکاک 
سطحی در دمای محیط المان‌ها دچار اعوجاج و به‌هم 
ریختگی شده و نتایج شبیه‌سازی را غیر قابل اطمینان 
می‌نمود. این موضوع در شکل (۵ - الف) نشان داده 
شده است. 

رابطه تتش کرنش از روابط فشار گرم در دماهای 
مورد استفاده در شبیه‌سازی و چندین نرخ کرنش 
استخراج شده است و به نرم‌افزار داده شده است. از 
آنجاکه نرم‌افزار اباکوس در مدول تعریف خواص خود 
مدلی از پیش تعریف شده همانند قانون جانسون - 
کوک يا زنر-هولومان ندارد و فقط بصورت نمودارهای 
تنش - کرنش وابسته به دما و نرخ کرنش داده می- 
پذیرد. داده‌های تست فشار بصورت تنش و کرنش 
حقیقی و وابسته به دما و نرخ کرنش وارد محیط نرم- 
افزار شدند. 

همانطور که مشاهده می‌شود. المان‌های سطحی 
دچار به‌هم ریختگی و اعوجاج شده‌اند. بنابراین جهعت 


سال بیست و هشتم شماره یک» ۱۳۹۵6 


رفع این مشکل. از تکنیک ۸3,۳ در شبیه‌سازی 02561 
بهره گرفته شد که نتایج آن در شکل (۵ -ب) نشان 
داده شده است. در این حالت مشکل به‌هم ریختگی و 
اعوجاج بر طرف شده است. 


جدول ۱ پارامترهای به‌کار رفته در شبیه‌سازی توسط 


نرم افزار آباکوس 
اتعات تین که کت ۱ 

ابعاد مقطع کانال که ۱۶۸۱۶ 

بیرعت ورن ۱ میلیمتر بر ثانیه 
ضریب اصطکاک ۰/۹۵ 

اندازه مش استفاده شده ۱ 

تعگالی نموثه /۲۷۰۰ 
مدول یانگ نمونه 2 1 
ضریب پواسون ۳۳ 


شکل ۵ شبیه سازی 1( الف) بدون استفاده از تکنیک ۰03 و 


ب) با استفاده از تکنیک ۸5 


نتایج شبیه‌سازی‌های بعدی نشان می‌دهد که استفاده از 
این تکنیک به‌منظور شبیه‌سازی فرایند فشردن در 
کانال‌های هم‌مقطع زاویه‌دار فقط در قالب ٩۰‏ درجه (یا 
کمتر) در دمای محیط در پاس‌های یک و بالاتر مورد 
نیاز است. در حالی که در موارد دیگر (دماهای بالاتر و 
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پاسهای کاری بیشتر)» بدلیل کرنش‌های اعمالی کمتر و 
همچنین کم شدن اثر اصطکاک سطحی در دماهای 
بالات استفاده از تکنیک "م۸ مورد نیاز نمی‌باشد. 

منظور از کاهش اثر اصطکاک با افزوده شدن دما. 
بهتر شدن روند سیلان ماده با فعال شده پدیده‌های 
ترمیم همانند بازیابی و تبلور مجدد دینامیکی است. 
بعبارت دیگر با افزایش دماء روند کارسختی بدلیل فعال 
شدن مکانیزمهای کارنرمی آهسته و متوقف می‌شود و 
این باعث سیلان بهتر ماده درون قالب با نیروی کمتر 
می‌گردد و این مساله پدیده اعوجاج مش‌بندی را مرتفع 
می‌سازد. 
شکل (1). تصاویری از توزیع کرنش نمونه در چهار 
زمان مختلف از شبیه‌سازی را نشان می‌دهد. با توجه به 
ابعاد نمونه و قالب و سرعت پرس. زمان مورد نیاز 
فرایند شبیه‌سازی» 1٩‏ ثانیه می‌باشد. 

همانطور که از شکل () مشخص است. توزیع 
کرنش از بالا به پایین نمونه کاهش می‌یابد که علت آن 
را می‌توان به اختلاف سرعت المان‌های سطوح بالا و 
پایین نمونه نسبت داد. به این معنی که المان‌های سطوح 
بالایی باید با سرعت بیشتری نسبت به المان‌های پایینی 
حرکت کنند تا هندسه تغییرشکل حفظ شود. در نتیجه 
میزان کرنش اعمالی در سطوح بالایی بیشتر می‌باشد. 


وید نکن( بقتان یتفن وگن نوی قوس 
برابر با ۱۳۹ می‌باشد. این مقدار با استفاده از رابطه 


تئوری ایواهاشی [26] به صورت زير محاسبه می‌شود: 


ات مس( کتک ]مد - 3 
7 | )مهم( )مد 1 
2 6 2 3 


بنابراین تطابق نسبتاً خوبی بین نتایج حاصل از 


شبیه‌سازی و روابط تشوری وجود دارد که نشان از 


صحت شبه‌سازی انجام شده می‌باشد. شایان ذکر 


تشرد مهندسی متالررژی و مراد 


شبیه سازی تغییرات ریزساعتا رآلومینیوم... 


است که تفاوت بین نتایج شبیه‌سازی و محاسبات 


تئوری به تاثیر ترکیبی شرایط اصطکاکی و خواص مواد 
شامل کرنش سختی و جهات کریستالوگرافی نسبت داده 
می‌شود که در محاسبات تئوری در نظر گرفته نمی‌شود 
[127 


۱ 


9 ات 
شکل توزیع کرنش در 02561 در الف) 205 ب) 2195 


پ) 246.25 و ت) 1-665 
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فاگ فا ان ان 


1 ۱ ۱ ۱ ۱ ۱ ۳ 
14 112 10 8 0 4 2 0 
(هت) و۵ ۳0۲ عع۵ه)کز۲ 


زب 
شکل ۷ توزیع کرنش در نمونه شبیه‌سازی شده. الف) موقعیت 
خط در قطعه ۵-5 ب) تغییرات کرنش در امتداد حط ۸-2 


تغییرات کرنش در راستای خط ۸ در مقطع میانی 
نمونه در شکل (۷) نشان داده شده است. همانطورکه از 
شکل (۷- الف) دیده می‌شود. تغییرشکل در راستای 
محور 7 تغییر چندانی نذارد و تقریباً یکنواخت می‌باشد. 
بنابراین» راستای ۸-3 برای نشان دادن توزیع کرنش بر 
حسب فاصله از سطح بالایی نمونه در این مقطع 
انتخاب شد. همانطور که در شکل (۷ - ب) مشاهده 
می‌شود. مقدار کرنش اعمالی روند کاهشی دارد. 
به‌گونه‌ای که توزیع کرنش در ۲/۶ میانی تقریبا 
یکنواخت بوده و در ۱/۶ بالایی و پایینی غیر یکنواخت 
می‌باشد. این اس همانطور که قبلا نیز اشاره شد. 
احتمالا بدلیل اختلاف میدان‌های سرعت در عرض 
نمونه است. شایان ذکر است که از منحنی‌های کرنش بر 
حسب جابجایی بعنوان داده ورودی به الگوریتم مونت 
کارلو به‌منظور پیش‌بینی متوسط اندازه دانه نهایی نمونه 
ان شده پس از انجام فرایند فشردن در کانال‌های 
هم‌مقطع زاویه‌دار استفاده شده است. 


سال بیست و هشتم. شماره یک» ۱۳۹۵6 


نتایج شبیه‌سازی تبلور مجدد به روش مونت کارلو 

ها رن ریا ان کجونه تقوم کی ون گنها 
هم‌سفطع زاویه‌دار پس از آنیل. به‌منظور شبیه‌سازی 
ریزساختار نمونه فشرده شده در کانال‌های هم‌مقطع 
زاویه‌دار. از کرنش‌های به‌دست آمده توسط شبیه‌سازی 
اجزاء محدود به عنوان ورودی برنامه مونت کارلو 
استفاده می‌شود, 

برای اعمال مسیر کرنش به الگوریتم مونت کارلی 
از تقریب تکه‌ایخطی آن استفاده می‌شود. به این 
منظور. همانطور که در قسمت قبل توضیح داده شد. 
سطح مقطع نمونه به ۱۰ منطقه تبدیل شده و برای هر 
منطقه. متوسط کرنش اعمالی محاسبه می‌شود. این 
کرنش‌ها در شکل (۸) به همراه نمودار کرنش اصلی 


نشان داده شده است. 


4 126 11,72 98 هو 7 و مه و2 م4 0 
(صص) ۱۵۵ 6۵ ععصها 1۳6 


شکل ۸ مسیر کرنش در سطح مقطع نمونه و تقریب 
تکه ای-خعطی آن 

تکامل ریزساختار پیش‌بینی شده طی تبلور مجدد 
در دمای ۱۲۳ کلوین پس از فرایند فشردن در کانال‌های 
هم‌مقطع زاویه‌دار با کانال ٩۰‏ درجه در دو پاس در 
شکل ٩(‏ - الف تا ح‌( نشان داده شده است. در این 
شکل‌ها؛ به‌منظور تمایز بین جوانه يا دانه‌های تبلور 
مجدد يافته و دانه‌های تبلور مجدد نیافته. دانه‌های تبلور 
مجدد يافته با رنگ سفید مشخص شده و مرزهای آن 
نیز با خطوط درشت‌تر نشان داده شده است. 

هبانظون. که مشاهده م‌شود جر مراسل اوه 
دانه‌های تبلور مجدد به صورت مجزا از یکدیگر رشد 
می‌کنند (شکل ٩‏ - الف تا پ). پس از آن, دانه‌های 
کار فوتق ور با نگل اه ی شرنه: ار ٩‏ 
(ت و ث)) و در نهایت. تبلور مجدد در کسر اندک 
نواحی تبلور مجدد نیافته رشد می‌کند (شکل ٩-ج).‏ 
پس از اينکه تبلور مجدد کامل شد. دانه‌های تبلور 


1 


مجدد يافته با سرعت کمی رشد می‌کنند (شکل 5ج 


و ج) 


۱۱65۳۸6: ۸6۵6۳23۱27 


(ث) (ج) 


شکل ٩‏ تکامل ریزساختار پیش‌بینی شده طی تبلور مجدد در دمای 


۳ کلوین پس از فرایند فشردن د رکانال‌های هم‌طع زاویه‌دار با 


کانال ٩۰‏ درجه در دو پاس 


همچنین مشاهده می‌شود که در برخی زمان‌ها. 
نواحی کوچکی از بخش تبلور مجدد نیافته توسط 
دانه‌های تبلور مجدد احاطه شده و اصطلاحا به دام 
می‌افتند (شکل ٩‏ (ث)). این پدیده با یافته‌های ایوانی, 


نشریه مهندسی متالویژی و مواد 


شبیه سازی تغییرات ریزساختا رآلومینیوم... 


که شبیه‌سازی مونت کارلو برای اکستروژن آلیاژ ۸1 
4.577-68 انجام داد. همخوانی دارد [28]. 

در شکل ( ۱۰ متوسط اندازه دانه ریزساختار 
نهایی پیش‌بینی شده در هر یک از مناطق ده‌گانه نشان 
داده شده است. همانطور که مشاهده می‌شود. متوسط 
اندازه دانه در مناطق ۱ تا ٩‏ تقریبا یکنواخت بوده در 
اندازه متوسط دانه رشد بیشتری داشته است. دلیل این 
کمتر » انرژی ذخیره شده در آن منطقه کمتر بوده و در 
نتیجه دانسیته جوانه‌زنی آن نیز کمتر خواهد بود که 
باعث می‌شود جوانه‌های تبلور مجدد ایجاد شده فرصت 
بیشتری برای رشد در نواحی تغییرشکل داشته باشند. 
این نتایج با [25] که در آن به بررسی ناهمگنی 
ریزساختار نهایی سیم تخت شده حاصل از نورد پس از 
آنیل پرداخته است. همخوانی دارد. با توجه به این 
نمودار می‌توان نتیجه گرفت که با افزایش ناهمگنی 
کرنش در سطح نمونه» ناهمگنی اندازه دانه‌ها در 
ریزساختار نهایی نیز بیشتر خواهد شد. بعنوان مثال, در 
صورتی که فرایند فشردن در کانال‌های هم‌مقطع 
زاویه‌دار در مسیر ۸ تا تعداد پاس‌های بالاتری ادامه 
یابده اختلاف کرنش بین سطح نمونه و پایین آن بیشتر 
شده و می‌توان انتظار داشتتت ناهمگنی ریزساختار نهایی 


نیز افزایش یابد. 


(صر) 5126 طنحتع عو۲2ه ۸ 


0عع؟۳ 


شکل ۱۰ متوسط اندازه دانه ریزساختار نهایی پیش‌بینی شده در هر 
یک از مناطق ده‌گانه 
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پیش‌بینی ریزساختار نمونه فشرده شده در کانال‌های 
هم‌سفطع زاویه‌دار در زوایای کانال مختلف . در این 
بخش» ریزساختار نهایی پیش‌بینی شده نمونه فشرده 
شده در کانال‌های هم‌مقطع زاویه‌دار در دو زاویه کانال 
۰ و ۱۲۰ درجه در دو پاس با یکدیگر مقایسه خواهد 
شد. در شکل (۰)۱۱ تقریب‌های تکه‌ایخطی کرنش بر 
حسب فاصله از سطحء که نتیجه بخش شبیه‌سازی اجزاء 
محدود می‌باشد» برای دو حالت ذکر شده نشان داده 
شلده است, هعانطور کهمشاهده هی شوده کرنش مربوط 
به زاویه کانال ۱۲۰ درجه کمتر می‌باشد. بنابراین» انتظار 
داریم ریزساختار نهایی آن دارای دانه‌های درشت‌تری 
باشد. زیرا برای این حالت انرژی ذخبره شده و در 
نتیجه دانسیته جوانه‌زنی کمتر خواهد بود. 

در شکل (۱۲» ریزساختارهای پیش‌بینی شده 


نشان داده شده است. متوسط اندازه دانه برای زاویه 


سوت ۱۲۰ات 


کانال ۰ و ۱۳۰ درجه به ترئیب برابر با ۱/۷ و 1۹۲ 
میکرومتر پیش‌بینی شده است. 

شایان ذکر است که در حالت زاویه کانال ۱۲۰ 
درجه به دلیل مقدار کرنش پایین. دانسیته جوانه‌ها در 
دو منطقه آخر (دورترین مناطق از سطح) صفر بوده که 
این عامل باعث می‌شود که ناهمگنی ریزساختار نهایی 
زیاد شود. در شکل (۱۳) متوسط اندازه دانه برای مناطق 
ده‌گانه نشان داده شده است. همانطور که از این نمودار 
نیز مشاهده می‌شود در حالت زاویه کانال ۱۲۰ درجه. 
ناهمگنی بیشتری مشاهده می‌شود. انحراف معیار 
داده‌های نرمال شده برای دو حالت زاویه کانال ٩۰‏ و 
۰ درجه به ترتیب برابر با ۰/۰۸ و ۰/۱۱ بدست 


می‌آید که نشان‌دهنده مطلب فوق می‌باشد. 
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(حصص) حرم) و ععصه)وزظ 


شکل ۱۱ تقریب‌های تکه ای-حطی کرنش بر حسب فاصله از سطح. برای فرایند فشردن د رکانال‌های 


هم‌مقطع زاویه‌دار در دو زاویه کانال ٩۰‏ و ۱۲۰ درجه در دو پاس 


«سیر۸۵64,9213 :۱۱۵65۵۵ 
(ب) 


شکل ۱۲ ریزساختار نهایی پیش‌بینی شده نمونه فشره شده د رکانال‌های همقطع زاویه‌دار در زاویه کانال الف) ٩۰‏ و 


ب) ۱۲۰ درجه در دو پاس 
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شکل ۱۳مومیط آندازه:دانه زیرساهتار این پیش نی ده در هر یک 


نتیجه گیری 


این پژوهش در دو بخش کلی شبیه‌سازی اجزاء محدود 

فرایند فشردن در کانال‌های هم‌مقطع زاویه‌دار با استفاده 

از نرم‌افزار آباکوس و پیاده‌سازی الگوریتم مونت کارلو 

در محیط برنامه نویسی متلب. به‌منظور پیش‌بینی 

ریزساختار نهایی نمونه پس از آنیل» انجام شد. از 

منحنی کرنش بر حسب فاصله از سطح نمونه بعنوان 

ورودی مونت کارلو استفاده شد. 

نتایج به‌دست آمده از شبیه‌سازی اجزاء محدود فرایند 

فشردن در کانال‌های هم‌مقطع زاویه‌دار بطور خلاصه 

بصورت زیر است: 

- نتایج شبیه‌سازی‌ها نشان می‌دهد که استفاده از 
تکنیک ۸ به‌منظور شبیه‌سازی اجزاء محدود 
فرایند فشردن در کانال‌های هم‌مقطع زاویه‌داره در 
قالب ٩۰‏ درجه (و يا کمتر)؛ در دمای محیط, در 
پاس‌های یک و بالاتر مورد نیاز است. در حالی که 
در دماهای بالاتر و پاسهای بیشت به دلیل 
کرنش‌های اعمالی کمتر و همچنین کم شدن اثر 
اصطکاک سطحی در دماهای بالات استفاده از 
تکنیک ۸ مورد نیاز نمی‌باشد. 

- نتایج توزیع کرنش نشان می‌دهد که کرنش از بالا به 
پایین نمونه فشرده شده در کانال‌های هم‌مقطع 
زاویه‌دار کاهش می‌یابد که دلیل آن, اختلاف 
سرعت المان‌های سطوح بالا و پایین نمونه 
می‌باشد. به این معنی که المان‌های سطوح بالایی 


شبیه سازی تغییرات ریزساختا رآلومینیوم... 
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ز مناطق ده‌گانه برای زاویه کانال ٩۰‏ و ۱۲۰ درجه 


باید با سرعت بیشتری نسبت به المان‌های پایینی 
همچنین نتایج برشتته. ۰۰ اف او یش ی 
ریزساختار با استفاده از الگوریتم مونت کارلو بطور 
خحلاصه به‌صورت زیر است: 
از کرنش‌های به‌دست آمده از شبیه‌سازی اجزاء 
محدود برای محاسبه توزیع انرژی ذخیره شده در 
شبکه اولیه ورودی به الگوریتم تبلورمجدد مونت 
کارلو استفاده شده و به اين ترتیب کوپل شبیه‌سازی 
اجزاء محدود و مونت کارلو انجام شد. 
با افزايش ناهمگنی کرنش در سطح نمونه, 
مشک دتانه داشها حن اسان هاش 
بیشتر خواهد شد. 

با افزايش زاویه کانال قالب فرایند فشردن در 
کانال‌های هم‌مقطع زاویه‌دان متوسط اندازه دانه 
نهایی پس از انیل افزایش می‌یابد. متوسط اندازه 
دانه برای زاویه کانال ٩۰‏ و ۱۲۰ درجه به ترتیب 
پرابر با ۲۷۶ و 2/۹۲ میکرومتر پیش‌بینی اشد: 
همچنین در حالت زاویه کانال ۱۲۰ درجه به دلیل 
مقدار کرنش پایین. دانسیته جوانه‌ها در دو منطقه 
آخر (دورترین مناطق از سطح), صفر بوده که اين 
عامل باعث می‌شود که ناهمگنی ریزساختار نهایی 


زیاد شود. 
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شماره بیست و هشتی ص. ۲۹-۳۷ ( ۱۳۸۳). 
.6۱۱۱ ووزااع۱۵0 راومه مس متام مناموممنم فصتالع۳۱۷۵0 ررژ.۲ ورمطمصت .9 
2000(۰) ,893-910 00۰ ر8 ۷۵۰ ر.ع۳ :50 
.2۵000-1ت ۲6۵۵ ۵ ممتملمصتی ‏ مانمو‌طجومن ‏ رظن جمقتعلصظ. ریویع) ‏ اقعاق) ‏ ونگیظا. عااما0و ...10 
1986(۰) ,1833-1845 ۲0۰ و9 ٩0۰‏ 34 ۷۵1 ,.۷/6۵۵۱۱ 2۸660 را 0تع 4صع صمتاهع ناه فیامع۲۱0۵۴9۵868 
0 ماگ0 طمه2 م۲6۵0 ۵۴ تصمصصمماع7ع0 مط؟ رب وفع ری عصطامیا ریک عفصتطه]۱۷ :11 
(1989) ,139-145 .۲0 م1 .0 ,37 ۷۵1 .۷6۵1۱ 2۸660 رخطمتاهااصتی عاناوجوی ود 24 توااهاعنهورم 
۱0-0 صا صمتاه2نماوو۵ع۲ عقاو ۵۶ حمتاقایصته جمتقصمانتة تقاتاای م1 ربا ۵ظ۲۱۵ ر.ل.0.۳ وعز1(27 .12 
(1999) ,1145-1150 .20 ,10 .۱0 ,40 ۷۵1 ,۷۵۱6۲۱۵10( ماوز 50 ,۸1050۳ 
,وه صنمتع ...۵ ومتاماحصته ماتام‌صطمن؟ ربق.ظ تصطهو وی فعتق) ر..نا 9۲۵10102 ۱۳ طمفعلصظ .13 
.(1984) ,783-791 .00 و5 ٩0۰‏ و32 ,۷۵۱ .۱۷۲۵ 2۸6۱۵ رخیناعظ۲ 
تعنامم صا منم صتقع ۵۶ مصصطاتمقاه مومتاقانصته ماجقن عظ۷]۵ ر.ظ 1۷۲۵۲۳۵۵60۷۵ ر.1۷ »۷۲۵۳۳۵ ,14 
2000(۰) 117-128 .00 رد ٩0.‏ وک ,۷۵1 .۲60۳۳10 ٩66.‏ ۲۷5 ر" اه 6۲] ۳2۵ 
1۱27۵1:00-1۰ ۲6۵0۵۵ ۵۶ جمتمااصته تمتناو‌جمن؟ نایم عمط ,۱۲.۲ جمعتع0ص ررو.ی) اوعتم ررژ. 2ازاما0ت5 .15 
.(1988) ,2115-2128 00۰ ,8 ۱۲0۰ ,36 .۷۵1 ,.۷/6۵۵/1 2۸660 رامع مجح صمتاجع1)ناه ونامعصعع ۲۱6۱6۲0 
مها تعطا هه ۲۵6۵0۷۵ رطمتاه2زالماووهعر من نی مصتاله۷۵0 ۳ ری۱۷ ۲رفصصجعاعط رک قع0و .16 
۰ بر ۱9۱866۲۱۱ ۵080 56۱6866 وام ۷/۵۵۵۲ رواهتتعاههه مرمع افتامعصع8۵مصمطصا مد مصعصعهع0 
2002(۰) ,153-160 0۳۰ ,332 
هام6۵ ۵۴ جمتاماته من .۱۷۲۳۲ وموتعصم ر.نایظ عاونا ر.ژینا 5101072 ,تایه 101166 .17 
٩0۰ 2, 00. 627-639, )1989(.‏ رد3 ,۷۵۱ .۱۷66۵11 2۸66۵" کاهاع 0عصملع0 «اصمکتصحجمهط ط1 
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۵۰ شبیه سازی تغییرات ریزساختا رآلومینیوم... 


103006606 ماگ ممتاج م۱6۵۵ ۵۴ همناع۱۷]۵0 ر.ط.۲ تعصر؟ ری م۱۷ ری رفجعتافط رک هومو .18 
(2001) ,2206 -2199 00۰ و32 ,۷۵۱ یه ۲۳۵۴۵۹۵۵606 عامزعع ۱۷۵ ۵۵ آمعا ۱۷/6۵۱9 رخ طمتاعصولع 

لقاع ۵۶ عصناهعصصه نس ممتانمی طونم مص همان فطع صمتاقلصنو؟ ری 0ق2عصنصم1۵2 .19 
(2008) ,309-312 00۰ ر43 ۷۵۱۰ ,66اه وا۱۵ ۷۵۵6۲ امصمتاماا۵۳۵ ,۲۵۸۲۳ زرط 2۲۵665560 

مامت احسا‌تتاومعنصه هه متاعصن! عمط مهصتلع۱۷]۵0 ر.] ماتقطعه2 ررظیعم) مهو ربظ صمصطوته ملق ,20 
۰ ,46 ,۷۵۱ ,.1۷/۵۵۵۲ ۸6۸۵ ر"طم۵1112۵01 ۲6۵۲/۵ ۵۴ ممتاهصدنه ملتهت موم -صما۵ امومع »اهاز فصتتنل 
1998(۰) ,4415-4433 .00 ,12 

7 ۱۳۵9۲ فرع 20 هصنمم0 2۲ ۷۵۲1 صتجتی م۱2۲۵ ور تمصع .۲۰۷ عنام ررل مصقهالزعو نی :21 
1980(۰) ,69-412 00۰ ,25 ۷۵۱۰ ,5616۳166 ۷۲۵۵6۲:۵5[ 

مصکوع۳ تعامج اعصصفد انوم عط ۵۶ هتم »تافقاظ؟ ررو.تا ۷۵۵۵ ریک لک ریک ۷ 98 ریل.۲۷۷ لک .22 
,360-368 .00 ,487 ۷۵۱ یگ ۱9۱66۲۱۵ ۵0۳0 5616866 کامز۱۵۱۵۲ ,"ماه صتصتصصاح 6061 2۲۵6۵9560 
2008(۰) 

۲۵۱0۲60 عازومموصم تاه اهاع 6061-له صه ۵۶ ۵۵۵۷10۲ ۲۵۵۵ ر ,)۰۲۰ 121982001 ور با .23 
.(1997) ,4797-4806 .00 رك1 .۲0 و45 ,۷۵۱ .۱۷۵66۲ ۸60 روعاهان‌تاتهم مصتصاه 

۵اه ,"۳۵06 صمتامتالموو۵ع۲ ماتقن) مامم۷ ۵۶ ممتای‌صننع۷ امتصمصنهص ۳۲ م1 ۲۷۷21266 .24 
2004(۰) ,262-267 00۰ 157-158 ۷۵۰ ,66۳۳00۵109 ودزعدع۴۳۵6 وام ۷۲۵۸۵۲[ 

عظ۷۵ 0صج عصعصصعام عاتصق عطا گم صمتاه12ا0؟ ره م۲1 امن رم اطع تصصتتی ,1۷۲ 20ظ2۳مدنجدم2ه .25 
۵و۵ ۵ امومع فنامعصمعمصممطهصا که ممتامعتالم یومع عمط عصتاماصته نم عم ماتهن) 
2007(۰) و765-773 00۰ و ,۷۵۱ ,ع06عه کا۲۱۵ع۷۵6 آم0 انامه 

مه اعصصعد-اهناوع ۵۶ ملوتمظ۳۳ظ ۰ ری.۰۲ 12282000 رب۱۷۲ مامصصعا ریت ۲۱۵۲162 ررگل ۷۷۵۵8۵ .۷ ۲۸۵۳۵50 .26 
.(143-146,)1996 .۵0 ر0ناه۱ ۷/۵6۵۲( 60۳۱26 رکاهاتعاههه 0عصتقتم مص-متاان ۵۶ مصتووعه۳0( عط ۶۵۲ مصنووع2۲ 

مه تعصری تام ۵۴ عععصمتالطا؟ رب تاباً رییط لاظ2 رب۷۲ لاأبا و.ل یا ونکلنظ ناعذا یآ لاو ون) لاب .ری 1618 .27 
و65 ۲قاتام‌صج . امصصعط. لفوع صا ممتتامتی ‏ معا 0ص صمتاقصمگعل . متافقام عطا صم رشت) 
2013(۰) ,1-9 00۰ ,6616066 کاه ۷۵66۲ آم0 نامه 

۶ حمتاهاصنه مه ماه 2نصععمصومط متس مامت ام وم منم ۵۶ مصتام۳۱۷۵0 ربنم ۳1۷28 .28 
۸۱-4.57,۳7-8 060ات)رع صز ریت بای تیم عفتقمم آهمطمتتمم ما مصنل۵ع1 ممتاهم1ه) ۲۱6۵6 عقاو امعتعصهط 
۶ انونه7 نا 1۳۵۶ مطا روممامصعع]. مهد معصمعتنی امتماع۱۷ ۵۶ صمماتهوم روتفعطا بلاط 2۱107۳ 
2010(۰) ,60۳۳001087 1[ 


تخر مهتلسی مالی زج و مراد سال بیست و هشتم شماره یک, ۱۳۹۵ 


